Nuclear Shadowing in Neutrino-Nucleus Deeply Inelastic Scattering 
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In the framework of the colhnear factorized pQCD approach we calculate the small-is process- 
dependent nuclear modification to the structure functions measured in neutrino-nucleus deeply in- 
elastic scattering. We include both heavy quark mass corrections (A/^/Q^) and resummed nuclear- 
enhanced dynamical power corrections in the quantity {^^"^ / Q^){A^^^ — 1) with evaluated to leading 
order in Qs . Our formalism predicts a measurable difference in the shadowing pattern of the struc- 
ture functions F^{xb,Q^) and -F3*(a;s,Q^) and a significant low- and moderate-Q'^ modification of 
the QCD sum rules. We also comment on the relevance of our results to the NuTeV extraction of 
sin^ 9w- 

PACS numbers: 12.38.Cy; 12.39.St; 24.85.-|-p; 25.30.-c 



I. INTRODUCTION 



Recent surprising results on sin^ 6w , reported by the 
NuTeV collaboration and based on a comparison of 
charged and neutral current neutrino interactions with 
an iron rich target , renewed our quest for understand- 
ing the nuclear dependence in neutrino-nucleus deeply 
inelastic scattering (DIS). A possibility that process- 
dependent nuclear shadowing might affect the NuTeV 
extraction of the Weinberg angle 9w was raised by Miller 
and Thomas 0- Although such scenario was considered 
unlikely by the collaboration Q, a systematic study and 
a clear understanding of the process-dependent nuclear 
effects in neutrino-nucleus scattering will strengthen the 
importance of the NuTeV result. 

Like all nuclear dependences in the physical cross sec- 
tions 4j, the small-XB shadowing in lepton- nucleus DIS 
has both process-dependent and process-independent 
contributions. While its universal part can be fac- 
torized in the leading twist nuclear parton distribu- 
tion functions (nPDFs), the DIS-specific modifications 
arise from the higher twist (or power) corrections to 
the structure functions 0, IE El- In this letter, we 
present a calculation of the process-dependent shadow- 
ing in neutrino-nucleus deeply inelastic scattering by re- 
summing heavy quark mass corrections, M^/(2p • q) = 
xbM^/Q^, and nuclear size enhanced dynamical power 
corrections, (CVQ2)(-^i/3 _ ^j^j^ ^2 ^ {f+^F\_), the 
gluon density in a large nucleus. The numerical value 
for the characteristic scale of higher twist £^ ex- 
tracted from DIS data on fi-A interactions is much 
less than in the region which is perturbatively accessi- 
ble. Therefore, we only evaluate to the leading order in 
as, while resumming the power corrections to all orders 
in (^V<3^)(^^^^ - !)■ Using = 0.09 - 0.12 GeV^ [|, 
our results provide a good description of the deviation 
between the Gross-Llewellyn Smith QCD sum rule 
adjusted for ©(as) scaling violations 9] and the existing 
data small Bjorken xb, the high twist compo- 

nents to the calculated structure functions F^{xb,Q^) 
and F^{xb,Q^) in neutrino-iron DIS qualitatively de- 
scribe the Iow-xb and low-Q^ suppression trend in the 



preliminary data, recently reported by the NuTeV col- 
laboration at DIS 2003 \u\ . 

In the next Section we briefly review the DIS kinemat- 
ics and coherence at small Bjorken j;^ . In Section III we 
demonstrate that at the tree level mass corrections and 
dynamical power corrections "commute" and their re- 
summation can be carried out in a closed form. We derive 
analytic expressions for the process-dependent nuclear 
modification to the transverse and longitudinal struc- 
ture functions in neutrino- nucleus DIS. In Section IV 
we predict the difference in the shadowing pattern of 
F^{xb,Q^) and F^{xb,Q'^), and give quantitative re- 
sults for the XB-, A- and Q^-dependence of the nuclear 
modification to the charged current v{v)-A DIS structure 
functions. We find sizable small- and moderate-Q^ cor- 
rections to the Gross-Llewellyn Smith QCD sum rule. In 
Section V we comment on the relevance of our results to 
the NuTeV extraction of sin^ 6w ■ Finally, we give our 
conclusions in Section VI. 



II. DIS KINEMATICS AND COHERENCE AT 
SMALL XB 

The charged current DIS cross section of a neutrino (or 
antineutrino) beam (fc) off a nuclear target {Pa)j as il- 
lustrated in Fig.QJa), probes three independent structure 

functions, F/^{xb, Q^) with i = 1, 2, 3 
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where the Bjorken variable xb = Q^/(2p • q) with 
p = Pa/ A, the exchanged W^-boson momentum q and 
its virtuality = —q^, and y = p ■ q/{p ■ k). In 
Eq. the "(— )" represents the sign for an antineu- 
trino beam, mjv — M^/A with nuclear mass M^, M\y 
is the VF-boson mass, and E is the beam energy. The 
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FIG. 1: (a) Illustration of the characteristic scales, A± — 
and Az'~' = 1/{xbP'^) in the boosted frame, probed by 
the virtual meson in DIS. (b) Multiple final state interactions 
of the struck quark with the partons from the nucleus at a 
fixed impact parameter b. 



often-referred longitudinal structure function, Fj 
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are predominantly interest in the small- a;^ region and 
neglect the target mass (rescaling) corrections 

The DIS cross section with an exchange of a W- or 
Z-boson of virtuality and energy i/ = Q"^ /{2mj\jXB) 
has an effective resolution in transverse area A± — 
which is much less than the nucleon size, and an uncer- 
tainty in longitudinal direction Az^~^ = 1/{xbP'^) with 
boosted nucleon momentum p+. If Az^^^ = 1/{xbP'^) > 
2ro (toat /p"*" ) or xb < xn = l/(2TOArro) 0.1, the neu- 
trino will coherently interact with more than one nucleon 
inside the nucleus, and probe the nuclear dependence at 
a perturbative scale H, |E Q ■ 



III. CALCULATING MASS AND DYNAMICAL 
POWER CORRECTIONS 



Elcctroweak charged and neutral current processes ne- 
cessitate a discussion of final state charm mass effects 
in neutrino-nucleus DIS even if the leading twist charm 
quark parton distribution is neglected, (j)c{x,Q^) = 
(j)c{x,Q'^) = Q. It is, therefore, critical to develop 
a systematic approach to the interplay of a heavy quark 
final state and the resummed nuclear enhanced power 
corrections discussed in 0. We define the boost invari- 
ant mass fraction 



XM 



Xb 



(2) 



and choose a frame such that p'' = p+n^ and = 
—XBP'^n'^ + Q"^ /{2xBP^)n'^ , where ri^ — [1,0, Oj_] and 
= [0,1, Oj_] specify the "-f" and "— " lightcone direc- 
tions, respectively. With a non- vanishing quark mass M, 
the Feynman rule for the final state cut line of quark 



momentum Xip -I- g in Fig. |2Ja) is 



Cut = 27r 



Xb 



■ p + M)S{xi - XB - xm) , (3) 



where 



xmP'^ + 



2xbP^ 



(4) 



with p^ — A'P . For M ^ we recover the known mass- 
less case, in which the scattered quark is moving along 
the "— " lightcone direction. A direct consequence of this 
Feynman rule is a tree-level coupling for longitudin ally 
polarized vector mesons oc 

AfVQ2 Contracting e^'' ^ 
with the charged current hadronic tensor W^y ^ yields: 



A 



F^ixB^Q' 



Af2 



U,D 



(5) 



-jFl^ixB.^ 



U.D 



\VUD?^<P^ {xB+XMo.Q^) 



D,U 



\'t^'I^d{=^b+xm^,Q') , (6) 



where the CKM matrix elements Vij parametrize the 
electroweak and mass eigenstate mixing with up-type 
quark U = {u,c,t) and down-type quark D = (d, s,6) 
,9], and (pf represent the flavor-i universal twist-2 par- 
ton distribution functions (PDFs) of a nucleon {A = 1) or 
a nucleus 0. Eqs. © and ® give a novel leading order 
(ofg) power suppressed (M^/Q^) quark mass contribu- 
tion to the ratio of longitudinal and transverse structure 
functions R{xb,Q^) ^ F^{xb, Q^)/Ff{xB,Q^) for both 
nucleons and nuclei. 

We calculate the nuclear enhanced dynamical power 
corrections in the lightcone A~^ — gauge. In this gauge, 
other than the initial-state contact-term contributions, 
all leading order nuclear enhanced power corrections are 
from final-state multiple gluon interactions of the scat- 
tered quark in a large nucleus shown in Fig. Efb) . To 
resum all order nuclear enhanced power corrections with 
a non-vanishing (anti)quark mass, we examine its prop- 
agator structure • For a quark momentum Xip + q 
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1 -p 



J -p + M 



{xb + xm) ± ie ' 
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where ±i, ±ie correspond to propagators to the left or 
right of the t = oo cut. In the Fourier space conjugate 
to Xip'^ the first term, free of Xi pole, is cx 5[y~). The 
operators in the hadronic matrix element that this con- 
tact term (—[—>) separates can be evaluated in the same 
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FIG. 2: (a) Tree level direct coupling of the exchange vector 
meson to the struck quark in charge current neutrino- 
nucleus DIS with massive charm final state, (b) Tree-level 
contribution to the nuclear-enhanced dynamical power cor- 
rections with heavy quark eflPects. 



nucleoli state Q. In contrast, the Fourier transform of 
the second term is oc 9{y^). Therefore, this pole term 
>) is the source of the yl^/"^ nuclear size enhancement 
to the power corrections. The operators that it connects 
in the multi-field multi-local hadronic matrix element can 
be long-distance separated and thus approximately eval- 
uated in different nucleon states . Alternative operator 
decompositions as well as other terms that arise from a 
formal operator product expansion (OPE) are sup- 
pressed by powers of the nuclear size. 

The case of massive final state quarks could be much 
more involved than the M — > limit. The complexity 
of the calculation stems from the potentially dangerous 
exponential growth of the number of terms coming from 
products of propagators, see Fig. ^h) . In our calculation 
we first observe that 'y-'p + M in the cut line, Eq. (O, and 
the numerator of the pole-term, Eq. Q, arise from an 
on-shell momentum p. Since the exchange gluons at the 
vertices connected by quark propagators are transversely 
polarized in the = physical gauge, Af^{y~)^^^ w 



' ■ p + M)j-^ ■ p + M) ■ ■ ■ oc -{p^ - M^) 7-' 







For the diagrams in Fig. |2Ib) only one alternating se- 
quence of short and long distance parts of the propaga- 
tors Eq. Q, initiated by the t — 00 cut, survives. There- 
fore, there must be an even number of gluon interactions 
between the cut line and any surviving pole term of a 
propagator in Fig.l^^b). We also note that 



■ 7 • p (7 • p ± M) 7 • p • • ■ oc 7 p^ 



\2xbP^ 



7 p+ 



leaves no mass dependence in the spinor trace of the di- 
agrams. The basic unit for two-gluon exchange with a 
net momentum fraction flow Xi — Xi-i and two-quark- 
propagators (one contact plus one pole) 0, see Fig.|2Kb), 
now reads: 



Unit = xb 



27r Xi — Xi^i — ie 

Xi^l-{xB+XM) + if- ' ^^^^ 

-^F^{^^) right 

Xi-(xB+XM)-il^ ' ° 



(8) 



In Eq. © the boost invariant Xi ~ P^y^ , the two cases 
correspond to a vertex to the left or right of the final state 
cut and F^{Xi) is given by the intra- nucleon two-gluon 
field strength correlator defined in [^: 



1 



dX, 



271 [p+y 



rF+"(A,)F/(AO0(A, -A,) . (9) 



•••7-p7 ^ ■ p ■ ■ ■ —p 7" 



= 0, 



We conclude that the dynamical nuclear-enhanced all 
twist contributions from the leading order in Feynman 
diagrams with a massive quark final state are identical to 
the massless case up to the substitution xb ^ xb + xm 
(rescaling) in the (5-function in the cut, Eq. ©, and the 
propagator poles, Eq. JSJ. Effectively, we have shown 
that the mass and nuclear enhanced power corrections 
"commute" and Xi = xb + xm for all "z" . This allows us 
to take all possible final state interaction diagrams and all 
possible cuts 0| to explicitly carry out the resummation 
of coherent high-twist contributions to neutrino-nucleus 
DIS structure functions, 



-F'(^{XB,Q^) ~ l^Dul^^i {xB + xm + XMu,Q^) ± E I^UdI^'I^O i^B + ^HT + XMs,Q^) , (10) 

^ \D,U u,D J 

^FliixB,Q^) « {2} I 5] \Vud\^cI,§ {xb + XHT + XMo,Q^) ± E l^Dul^^D {^B + ^HT 
\U.D dXJ 



+ XM^,Q^) ■ (11) 



In Eqs. (|l()(l and Ullfl the "±" signs refer to Fi (parity conserving) and F3 (parity violating) transverse struc- 
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ture functions, respectively. The factor "{2}" gives the 
standard normahzation for Ft, only 9] and the isospin 
average in the PDFs over the protons and neutrons in 
the nucleus is impHcit. In Eqs. ifTUIl and ((TT|l xht is 
the momentum fraction shift (rescaling) induced by nu- 
clear enhanced dynamical power corrections and derived 
in Ref. 6]: 



= (12) 



where f ^ represents the effective scale for the dynamical 
power corrections. To the leading order in ag it is given 
by 



? = ^72 {P\P \P) 



8r, 



(13) 







where {p\F'^\p) depends on the small- a; limit of the 
gluon distribution in the nucleon/nucleus 3- While 



xn ~ 0.1 is the limiting value for the onset of coher- 
ence, at XA ~ 1/(2tojv?'o^^^'^) < the exchange vector 
meson already probes the full nuclear size, see Fig.^ To 
first approximation, the function 



, xb > xn 

, xa<xb <XN 

1 , Xb < XA 



(14) 



in Eq. 112|l represents the interpolation between the two 
regimes based on the uncertainty principle [2(^. The nu- 
clear enhancement factor {A^^^^ — 1) in Eq. (|12|l comes 
from the integration J dXi in Eq. (jSJ, the lower limit of 
which was chosen such that the effect vanishes for the 
proton {A = 1) case. 

Including the dynamical power corrections, the longi- 
tudinal structure functions in Eqs. (O and © become: 



-Fl^ixB,Q') 



F^^^\xB,Q^) + J2\Vnu\ 

D,U 



M + ^(2 



Ml 



'^D {^B + XHT + XMu ! Q^) 



U,D 

(LT), 



Ml £2 
— + — ( 2 - 

n2 ^ n2 ' 



Ml 



Q^ + Ml 



bfj {xB + a;HT + a; Mo . Q^) 



(15) 



:Fl^{xB,Q^) « Ft^\x,Q^)+Y,\VuD 



U,D 



Ml , £2 ^^2 
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Q2 ' Q2 Q2 + 



D 



{xB + XUT + XMo I 



D,U 



^4 , e 



2 - 



M}, 



{xB + XWI + XMfj , Q^) 



(16) 



In Eqs. H15|l and (|16|1 we include the 0{as) leading twist 
longitudinal structure functions f'^"^\x,Q'^) since 
they are of the same order as the leading power. For 
numerical evaluation in the next Section we consider two 
quark generations, U — {u,c) and D = (d, s), and use 
IKdP = \Vcs\^ = cosHc = 0.95, Wus? = \Vcd? = 
SIT? 6c = 0.05 with Cabibbo angle 9c 19]. The u, d and s 
quarks are treated as massless and the charm quark mass 
is set to Mc = 1.35 GeV Q. 



IV. HIGH TWISTS, SHADOWING AND QCD 
SUM RULE 

We first quantify analytically the differences in the 
"shadowing" pattern induced by valance and sea quarks, 
neglecting the charm mass effects that are shown to be 
small below. For isoscalar-corrected {Z = N = A/2) 
target nuclei we average over neutrino- and antineutrino- 



initiated charged current interactions, F/^{xb,Q^) — 
{F^^{xB,Q'^) + F^'^{xB,Q'^))/2. In the leading-order 
and leading twist parton model F^[xb, Q^) measures the 
valance quark number density with (pvaiXx) oc 
at small x. F^{xb,Q^), a singlet distribution, is pro- 
portional to the momentum density of all interacting 
quark constituents and for xb <SC 0.1 is dominated 
by the sea contribution, 4>seaix) oc x~°''"^. There- 
fore, the a;B-dependent shift from dynamical nuclear en- 
hanced power corrections, Xht in Eq. itT^ . generates 
different modification to F^{xb,Q^) and F^{xb,Q^)- 
Let i?"^^"^ , (xr, Q^) be the shadowing ratio determined 
from nuclei A and A' (for example ^^Fe to ^D) in Fj and 
F3 . If the scale of high twist corrections <C [a| and 
Xb < Tahi{xA, XA'), we expand the PDFs to first order 
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0.0001 



FIG. 3: The predicted nuclear modification for isoscalar- 
corrected ^^C, ^^Fe and ^°*P6 to the neutrino-nucleus DIS 
stricture functions F^(xb, Q^) (top) and xbF^{3:b, Q^) (bot- 
tom) versus Bjorken xb (left) and (right). The bands 
correspond to = 0.09 - 0.12 GeV^ 



in xht to obtain: 



(17) 



Since a^ai. ~ 0.5 and asea ~ 1 vary slowly with Q^, 
Eq. (|17|) predicts a measurable difference in the nuclear 
shadowing for the structure function F2 (Fi) in compar- 
ison to F3. 

Figure O shows the modification to the DIS structure 
functions from Eqs. iPH-QSl) for large nuclei C^C, ^^Fe 
and ^°^P6) relative to the deuteron, calculated with the 
CTEQ6L parton distribution functions j23|. The bands 
correspond to a scale for power corrections — 0.09 — 
0.12 GeV^ extracted from the analysis ||| of the NMC 
and E665 data 0. The transition 0.01 < xb < 0.1 region 
represents the onset of coherent interactions, Eq. (|14|l . 
and the modest a;B-dependence for xb l£ 0.01 is driven by 
the change in the local slope of the PDFs. The right pan- 
els show the dependence of the modification from the 
resummed nuclear-enhanced power corrections, which is 
noticeably stronger than the DGLAP evolution of lead- 
ing twist shadowing in the nPDFs The difference 
in the suppression pattern of F^ and XbF^ in Fig. is 
qualitatively described by Eq. p7ll . In contrast, it has 
been suggested in the framework of a Gluber-Gribov ap- 
proach [23 that the suppression of the non-singlet distri- 
bution may be significantly larger than of the singlet one 
(■^fei^ai > 1 for A > A'). Such distinctly different pre- 
dictions should be testable in the future i/-Factory exper- 



iments, for example at the Fermilab NuMI facility |24| . 

Although the current v(i>)-A DIS measurements are 
mostly on nuclear targets [23, these data lack the nec- 
essary atomic weight systematics to identify small- ccb 
nuclear shadowing. We do, however, note that our re- 
sults provide a consistent explanation of the observed 
small- and moderate-Q^ power law deviation at small- 
xb of the preliminary NuTeV data on F^{xb,Q^) 
and xbF^{xb, Q^) from the next-to-leading order lead- 
ing twist QGD predictions using MRST parton distribu- 
tion functions |26j| . 

The latest global QGD fits include v{h')-A DIS data 
without nuclear correction other than isospin 21]. Such 
analysis would tend to artificially eliminate most of the 
higher twist contributions discussed here due to a trade 
off between the power corrections in a limited range of 
and the shape of the fitted input distributions at Qqj 
especially within the error bars of current data. An ef- 
fective way to verify the importance of the nuclear en- 
hanced power corrections for neutrino-nucleus deeply in- 
elastic scattering is via the QGD sum rules, in particular, 
the Gross-Llewellyn Smith (GLS) sum rule Q 



GLS 



dxf 



1 

2^ 



^bF^^) 



(18) 



At tree level Eq. (|18|l counts the number of valance 
quarks in a nucleon, 5gls = 3. Since valence quark num- 
ber conservation is enforced in the extraction of twist-2 
nucleon/nucleus PDFs, the adjustments of input parton 
distributions can alter their shape but not the numerical 
contributions to the GLS sum rule. 

The effect of scaling violations can modify ^glSi £^rid 
at 0{as) [3 



Agls = ^ (3 



a.(Q2) 



(19) 



Loop contributions to the GLS sum rule are known to 
0{al) 01 . Although power corrections can also modify 
the shape of nucleon structure functions, recent precision 
DIS data on both hydrogen and deuterium targets from 
JLab indicate that effects from higher twist to the 
lower moments of structure functions are very small at 
as low as 0.5 GeV^, which confirms the Bloom-Gilman 
duality [l^. A recent phenomenological study [s^l also 
suggests that power corrections to the proton F2{xb,Q'^) 
have different sign in the small- and large-x^ regions and 
largely cancel in the QGD sum rules. 

On the other hand, the coherence between different nu- 
cleons inside a large nucleus is only relevant for a;^ < xn. 
The suppression of structure functions at small Bjorken 
Xb in Fig. |3| caused by the nuclear enhanced dynamical 
power corrections, cannot be canceled in the moments 
and further reduces the numerical value of S'gls- Fig- 
ure ^ shows a calculation of Agls from Eqs. (|10|) and 
(|ll|l for ^^Fe. While the effect of charm mass is seen to 
be small relative to as/ir, for ^ 1 GeV^ nuclear en- 
hanced higher twists may contribute as much as ~ 10% 
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FIG. 4: Top panel: Agls calculated to 0{as) with charm 
mass (Mc = 1.35 GeV) effects (stripes) and Mc + resummed 
power corrections (band). Data is from CCFR, CHARM and 
IHEP-JINR 



to Agls ■ Their behavior is consistent with the trend 
in the current data jT(il |. For comparison, we found no 
deviations for this sum rule induced by the EKS98 scale 
dependent parameterization of nuclear effects '22'| , which 
is again a consequence of the valance quark number con- 
servation in leading twist shadowing. 



IMPLICATIONS FOR EXTRACTION OF 



Based on a comparison of charged and neutral cur- 
rent neutrino interactions (separated on the basis of event 
topology) with an iron-rich heavy target, the NuTeV col- 
laboration reported a measurement of 



q(on- 



■shell) 



0.2277 ± 0.0013(stat.) ± 0.0009(syst. 



neglecting the very small top quark and Higgs mass cor- 
rections. This result is approximately 3 standard devi- 
ations Pi above the Standard Model (SM) expectation 
value sin^ 9w = 0.2227 ± 0.0004. The NuTeV's result 
was derived from a quantity that is a close approxima- 
tion to the Paschos-Wolfenstein relationship 



,.NC 



R- = 



(^) 



(^) 



rCC 



(^) 



1 2 

sin 6w ■ 

2 



(20) 



Corrections to Eq. (|20|l include higher order and nonper- 
turbative QCD effects, higher order electroweak effects 
and nuclear effects. Based on a QCD glob al analysis of 
parton structure, Kretzer et al. argue [3 If that the un- 
certainties in the theory which relates R~ to sin^ 6\y are 
substantial on the scale of the precision NuTeV data and 
suggest that the sin^ 6w measurement, NuTeV dimuon 



data and other global data sets used in QCD parton 
structure analysis can all be consistent within the SM. 

Because a heavy target was used, several nuclear ef- 
fects can enter the cross sections to influence the extrac- 
tion of sin^ 9w |32 . Since nuclear enhanced power cor- 
rections were not included in NuTeV's analysis, Miller 
and Thomas pointed out that nuclear shadowing from a 
vector meson dominance (VMD) model could affect the 
charged and neutral current neutrino scattering differ- 
ently, and therefore change the predictions for the ra- 
tios of neutral current (NC) over charged current (CC) 
cross sections, R"'^^^ = cr^^(z^(p))/cr^^(i/(z?)), and the 
extraction of sin^ 9w ^ ■ The NuTeV collaboration ar- 
gued 0, 113 that such possibility was considered unlikely 
because R~ has little sensitivity to process-dependent 
nuclear effects. 

In this letter we calculated the process-dependent nu- 
clear effects in the neutrino-nucleus differential cross sec- 
tions, Eq. C^l, in the perturbatively accessible DIS region. 
Our predictions on the nuclear modification to the v{v)- 
A structure functions in Fig. O should be relevant for 
between 1 and 10 GeV^. While for the mean {Q'^)y — 
25.6 GeV^ and (Q^)^ = 15.4 GeV^ the effect of dynam- 
ical power corrections is small, a large fraction of the fi- 
nal data sample cover the xb < 0.1, < 10 GeV^ range 
where shadowing can be as large as ~ 20%. We note that 
the NuTeV measurement constitutes a ~ 2% increase in 
the value of sin^ 9w relative to the SM, or equivalently 
^ 4% reduction of the expected total neutrino-nucleus 
cross section. Including shadowing into the expected to- 
tal cross sections will certainly reduce the discrepancy 
of sm^9w. However, without knowing the nuclear en- 
hanced power corrections to the structure functions at 
< 1 GeV^, and the detailed Monte Carlo simulation 
of event distributions, it is difficult to estimate the pre- 
cise corrections to the extraction of sin^^vi^. We, nev- 
ertheless, note that at small Bjorken xb, the calculated 
nuclear structure functions F^{xb,Q^) and F^{xb,Q^) 
in neutrino-iron DIS qualitatively describe the Iow-Xb 
and low-Q^ suppression trend in the preliminary data, 
presented by the NuTeV collaboration at DIS 2003 [l^. 



VI. CONCLUSIONS 

In the framework of the perturbative QCD coUinear 
factorization approach [sj l34l |. we computed and re- 
sumed the tree level perturbative expansion of nuclear 
enhanced power corrections to the structure functions 
measured in inclusive (anti)neutrino-nucleus deeply in- 
elastic scattering. We demonstrated that these cor- 
rections commute with the final state heavy quark ef- 
fects and identified the new contributions to the lon- 
gitudinal structure function F^{xb,Q^)- Our calcu- 
lated Q^-dependent modification to the Gross-Llewellyn 
Smith sum rule agrees well with the existing measure- 
ments on an iron target [Tol |. Our approach predicts a 
non-negligible difference in the small- a:^ shadowing of 
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the structure functions F^{xb,Q^) {F^{xb,Q^)) and 
F^{xb,Q^), which is consistent with the trend in the 
prehminary NuTeV data^J- Although our results, valid 
in the perturbative region, are unlikely to have an imme- 
diate impact on the NuTeV's extraction of sin^ Ow^ the 
predicted xb-, Q^-, and A-dependence of the structure 
functions in the shadowing re gion can be tested at the 
future Fermilab NuMI facility |2J| . 
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